I 


A • ^ 


NASA TECHNICAL NOTE 


>o 


i 




LOAN COPY: RETU 

AFWL (DOOLSS 
KIRTLAND AFB, 





ANALYTICAL AND EXPERIMENTAL 
INVESTIGATION OF THE COAXIAL 
PLASMA GUN FOR USE 
AS A PARTICLE ACCELERATOR 

/ 

I 

I 

by Edward L. Shriver ' 

George C. Marshall Space Flight Center 
Marshall Space Flight Center , Ala. 35812 



i 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. 


APRIL 1972 


TECH LIBRARY KAFB, NM 


TECH LIBRARY KAFB, NM 



1. REPORT NO. 2. GOVERNMENT ACCESSION NO. 

NASA TN D-6687 

4. TITLE AND SUBTITLE 

ANALYTICAL AND EXPERIMENTAL INVESTIGATION OF THE 
COAXIAL PLASMA GUN FOR USE AS A PARTICLE 
ACCELERATOR 

7. AUTHOR (S) 

Edward L. Shriver 

3. RECIPIENT’S CATALOG NO. 

5. REPORT DATE 

April 1972 

6. PERFORMING ORGANIZATION CODE 

8. PERFORMING ORGANIZATION REPORT » 

9. PERFORMING ORGANIZATION NAME AND ADDRESS 

George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 

10. WORK UNIT NO. 

11. CONTRACT OR GRANT NO. 

13. TYPE OF REPORT & PERIOD COVERED 

Technical Note 

12. SPONSORING AGENCY NAME AND ADDRESS 

National Aeronautics and Space Administration 
Washington, D. C. 20546 

14. SPONSORING AGENCY CODE 


15. SUPPLEMENTARY NOTES 


16. ABSTRACT 

The coaxial plasma accelerator for use as a projectile accelerator is discussed. The accelera- 
tor is described physically and analytically by solution of the circuit equations, and by solving 
for the magnetic pressures which are formed by the j xB forces on the plasma. 

It is shown that the plasma density must be increased if the accelerator is to be used as a pro- 
jectile accelerator. Three different approaches to increasing the plasma density are discussed. 

When a magnetic field containment scheme was used to increase the plasma density, glass beads 
of 0.66 millimeter diameter were accelerated to 7 to 8 kilometers per second velocities. Glass 
beads of smaller diameter were accelerated to more than twice this velocity. 


17. KEY WORDS 

Particle Acceleration 
Plasma Density 
Glass Beads 


19. SECURITY CLASS1F. (of this report) 


20 . 


Unclassified 


16. DISTRIBUTION STATEMENT 


SECURITY CLASS1F. (of this page) 


21 . 


NO. OF PAGES 


22 . 


Unclassified 


77 


PRICE 


$ 3.00 


For sale by the National Technical Information Service, Springfield, Virginia 22151 




TABLE OF CONTENTS 


Page 

CHAPTER I. DESCRIPTION OF THE INVESTIGATION 1 

Introduction i 

Description of Experiment i 

Circuit Constant Determination 5 

CHAPTER H. THEORY OF THE COAXIAL PLASMA 

ACCELERATOR 8 

Development of the Circuit Equations 8 

Solution of the Circuit Equations 13 

CHAPTER HI. THE LOW-FREQUENCY MAGNETIC FIELDS IN A 

DISC ARISING FROM RADIAL CURRENTS 15 

Development of the Expression for the Magnetic 

Field 15 

The Magnetic and Thermal Pressures on the 

Disc 27 

CHAPTER IV. DISCREPANCIES BETWEEN THEORETICAL 

CALCULATIONS AND EXPERIMENTAL. RESULTS . . 34 

CHAPTER V. THE ENERGY AND MOMENTUM OF THE PLASMA . . 37 

The Energy Contained in the Plasma 37 

Transfer of the Plasma Momentum to a 

Projectile 42 

CHAPTER VI. SUMMARY AND CONCLUSIONS 45 

Summary 45 

Conclusions 46 

APPENDIX A. INDUCTANCE AND CAPACITANCE OF COAXIAL 

CONDUCTORS 47 

APPENDIX B. LCR CIRCUIT THEORY 52 



TABLE OF CONTENTS (Concluded) 


Page 


APPENDIX C. COMPUTER PROGRAM 57 

APPENDIX D. BEHAVIOR OF AN IONIZED DISC CONDUCTOR IN 

A COAXIAL PLASMA ACCELERATOR 67 

REFERENCES 70 

BIBLIOGRAPHY 72 



I 


LIST OF ILLUSTRATIONS 

Figure Title Page 

1. System Schematic 2 

2. Ringing Frequency Current Versus Time 6 

3. Schematic Equivalent Circuit 9 

4. Resistance Element of a Disc Conductor 13 

5. Cylindrical Coordinates of a Circular Cylindrical Disc ... 17 

6. Plasma Plume Leaving the Muzzle of the Coaxial 

Accelerator 32 

7. Coaxial Plasma Generator 33 

8. Theoretical Current Versus Time Curve 35 

9. Current Versus Time When Firing Coaxial Gun 35 

10. Distance Versus Time 38 

A-i. Hollow Conducting Cylinder 47 

C-l. Computer Main Program 59 

C-2. Computer SETONE Program 60 

C-3. Computer SIMPSO Program 61 

C-4. Computer SETTWO Program 62 

v 



LIST OF TABLES 


TABLE 

1 


TITLE 

Physical Characteristics of Equipment 
in Coaxial Accelerator Circuit 


PAGE 

4 



CHAPTER I 


DESCRIPTION OF THE INVESTIGATION 


Introduction , , . t!i . '[ 

Projectiles with velocities up to 10 kilometers/second are produced in 
the laboratory [1] and are used to study the physical processes which take 
place in materials impacted with very high speed projectiles. It. would also 
be of interest to produce projectiles with velocities of 20 to 30 kilometers/ 
second to simulate meteorite impacts. Such high velocity projectiles cannot 
be produced in a light gas gun like the one used in production of the 
10-kilometer/second projectiles. 

It has been suggested that a high velocity plasma generator might be 
used to accelerate projectiles to the 20-kilometer/second velocity range; this 
method will be discussed. 


Description of Experiment 

The apparatus for this experiment consisted of a capacitor bank stor- 
ing 34. 56 kilojoules of electrical energy at 16 kilovolts, a set of coaxial leads, 
an ignitron-switching bank, a coaxial plasma generator, a tubular-vacuum 
range, and associated instrumentation required to operate this system 
(Fig. 1). All of the components consisted of commercially available items 






with the exception of the coaxial plasma gun which was made of materials 
readily available commercially and assembled in the laboratory. The 
dimensions and physical characteristics of the system components are listed 
in Table 1. 

The coaxial gun was loaded by placing an aluminum foil disc, 

0. 053975 meter in diameter and 6. 35 x 10 -6 meters thick, in the barrel of the 
gun, 16. 5 centimeters from the muzzle, so that electrical contact existed 
between the foil and both the center and outer electrodes of the gun. The gun 
was mounted in one end of the tubular vacuum range and the diagnostic instru- 
mentation was mounted on the vacuum range. 

The capacitor bank was charged from a console which provided a 
rectified high voltage power supply and a safety -dump circuit which could, 
if necessary, bleed the charge off the capacitor bank. The capacitor bank 
was discharged through the coaxial plasma gun by triggering the ignitron bank 
through a thyratron-actuated trigger circuit. 

When the capacitor bank was discharged, an emf was impressed 
across the foil disc causing it to conduct current, evaporate, and ionize by 
resistive heating of the conducting foil and colliding of the aluminum atoms. 
Since there was a magnetic field-around the center electrode and since the 
current in the foil was perpendicular to this field, the charged particles were 
driven down the barrel of the plasma gun by the j x B forces. This mecha- 
nism produces a very high velocity gas which may be used to accelerate 
projectiles. 
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Circuit Constan t Determination 

The capacitance, inductance, and resistance of the circuit must be 
known before any calculation can be made. Since the capacitance of the bank 
was known and was orders of magnitude higher than any of the capacitances of 
either the coaxial cables or the gun, which were in parallel with the bank, the 
capacitance of the circuit was considered to be that of the capacitor bank. 

This was not the case for the circuit inductance, however, since the 
capacitor bank was very low in inductance compared to the other elements of 
the circuit. The inductance contributed by circuit components of coaxial con- 
figuration may be calculated (Appendix A) , but the inductance of the ignitrons 
and the terminal connections must be calculated from ring frequency measure- 
ments of the circuit (Appendix B and Fig. 2) where 


R 


f = 

o 27 r 


L C 
o 


4L 2 
o 


= 13 172 hertz 


Also, by measuring the current amplitude in the ringing circuit, the circuit 
resistance may be calculated from the expression (Appendix B) : 



The current curve was obtained by placing a test coil of a few turns of 
wire near the leads from the ignitron bank. [No attempt was made to calibrate 
the current response since a known current of sufficient strength for 
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FIGURE 2. RINGING FREQUENCY CURRENT VERSUS TIME 
calibration was not available. However, the frequency and relative amplitudes 
of the current peaks can be obtained from the current curve (Fig. 2).] A 
short circuit was introducted across the plasma gun, and a voltage of 3000 
volts was impressed across the capacitor bank. The bank was then discharged 
and the circuit underwent a damped oscillating discharge. 

The ringing frequency from this curve was determined (Appendix B) 
by measuring the time interval between zero crossings on the curve. This 
gives the half period: T = — . 

O Z I 

o 

Thus, the frequency, capacitance, inductance, and resistance of the 
circuit were determined: 

f = 13 172 hertz (Ringing frequency without accelerator) 
o 

C = 270 microfarads (Circuit capacitance without accelerator) 

L = 0. 4157 microhenry (Circuit inductance without accelerator) 
o 

R = 6. 881 x 10" 3 ohms (Circuit resistance without accelerator) 
o 
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Since the plasma gun was only a section of coaxial transmission line 
for the current-carrying portion, the inductance may be calculated: 


A = In henry/meter 


A. = 2 x 10" 7 In 


3. 8 


0. 635 

A = 0. 3578 microhenry /meter ( accelerator inductance) 

Also, the resistance of the gun was measured and was found to be 0. 003 ohm. 
The gun was 0. 67 meter long, giving a value of Rj = 0. 00448 ohm/meter. 


7 


CHAPTER H 


THEORY OF THE COAXIAL PLASMA ACCELERATOR 


Development of the Circuit Equations 
When the aluminum disc was inserted into the gun and the capacitor 
bank was discharged through the circuit, the discharge was not simply a 
damped periodic oscillation. The motion of the plasma down the gun barrel 
caused the current-carrying length of the gun to vary, thereby changing the 
circuit inductance and resistance. 

The equivalent circuit is shown in Figure 3. It can be seen that the 

circuit consists of fixed capacitance C, fixed inductance L , and fixed 

o 

resistance R^, all in series with a varying inductance Ax and a varying 
resistance Rjx. 

The potential drops around the circuit are given by Kirchoff's law [2] : 


-$7 [ (L + Ax)I] + (R +R lX )I + V^-V = 0 , (1) 

dt o o T 

where V^, is the potential drop across the plasma caused by the energy loss 

d 

which is used to heat the plasma, — £ ( L + Ax) I] represents the potential 

at o 

across the inductors, and V is the potential across the capacitor. 


V 




I dt 


( 2 ) 
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FIXED INDUCTANCE VARIABLE INDUCTANCE 



FIGURE 3. SCHEMATIC EQUIVALENT CIRCUIT 


After differentiation, equation (1) becomes 

(L 0 + Xx) (lf) + [ R 0 + R ' x + X I s -] i + v t - V = 0 . (3) 

The energy balance may be expressed as 

T CV o = + l‘ L „ + Xx > 12 + / <K C -R.x> i^dt 


+ / v Idt+ jcv 2 
0 z 


(4) 


t 

The term J V I dt is the time integral of the power absorbed in heating the 
0 

plasma, and the kinetic energy of the moving plasma is 
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Differentiating equation (4) with respect to time produces the expression 


0 = in 




. + * V*Cv(f) . 

Dividing this expression by I produces a new equation: 

* = f(i!) ♦ < L o (-f ) + Mfr) 

+ (Ro + Rix)i + Vt+ ^v (^) , 

dV I 

but from equatio^ ( 2) , = - — and 

m /d 2 x\/dx\ , . / d I \ I / dx \ 

0 = t Ur] Ur ,L o tW U) t i 1 ra 


+ (R o + R!x)I + V t -V . 

Now, substract equation (3) from equation (5) to obtain the following 
equation: 


(5) 


f (£)(£) 


or 


d^ 

dt 2 


I 2 


2m 


( 6 ) 


If the circuit constants and initial conditions are known, the coupled 
set of equations ( 3) and ( 6) may be solved simultaneously to obtain the 
potential drop across the capacitor, the plasma velocity, and the current. 
The two equations are 
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and 


d z x 

dt 2 


I 2 

2m 


(L + Xx) 
o 


f + ' [ E o + R '* + *(f)] 


I + V T -V=0 . 


These equations, without the resistance terms or plasma heating term, have 
been solved by others [3,4,5] using dimensionless variables. The dimen- 
sionless set of equations was 

^ - ki* 
dT 2 


and 


1-fidT = i + £ — , 

J dr s dT 

(s is an equivalent distance) and 


C 2 V 2 X 


where the constant k = 


2ms 


L + Xx = Xs (x is the length of the gun barrel on the capacitor side of the 
o o o 

aluminum disc). It was found that a value of k = 1 most nearly fits the 
observed discharge conditions. Using this value for k and solving for m, 
it was found that the mass of the accelerated plasma is 


m = 


C 2 V 2 X 2 
o 

2(L + Xx ) 
o o 


(7) 


Substituting this mass expression into equation (6) , the differential equation 
becomes 


d 2 x 

dt 2 


L + Xx 

o O 1 J 2 


c 2 v 2 X 
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To solve the coupled set of equations , the voltage arising from the 
resistance of the plasma, V^, , must be expressed as an analytic function of 
the current. 

If it is assumed that the resistivity of aluminum remains constant 
during the heating and ionization process , the voltage across the foil can be 
calculated. The resistivity of the aluminum at the melting point is [ 6] : 
p = 10. 7 x 10 -6 ohm - cm 

The resistivity is a nearly linear function of the temperature throughout the 
solid phase, but no information is available for the liquid phase. In any case, 
the resistivity will probably not change by orders of magnitude as long as the 
density remains constant. Then, the resistance of the aluminum disc can be 
calculated where the resistance is 



l is the conductor length and A is the cross-sectional area 
For a conductor of annular cross section of length 1, 
and radius a, the resistance is (Fig. 4): 
da 


dR 2 — p 


27rla 


of the conductor, 
thickness da, 


and 


r 2 



a 

max 


a . 
min 


Then, for the disc used, R 2 = 4. 798 x 10 3 ohms. 

The voltage drop across the disc will then be 

v T = I R 2 . 
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Solution of the Cir cuit Equations 

These equations were solved on the Sperry-Rand 1108 computer using 


a fourth-order Runge-Kutta program. When solving the equations, the 
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simultaneous solution was used until the value of x was equal to the length 


of the gun barrel. At this time, only the damped 
applicable, since ^ is constant and = 0. 

conditions was used (at time t = 0) (Appendix C) 

x = x = 0. 505 meter 
o 

V = V = 1. 6 x 10 4 volts 
o 

L = 4. 157 x 10 -7 henry s 
o 

X = 3. 578 x 10 -7 henrys/meter 
R = 6. 881 x 10 -3 ohms 


o 

Rl = 

R 2 = 

C = 

dx 

dt 


4. 48 x 10 -3 ohms 
4. 798 x 10 -3 ohms 
2. 7 x 10~ 4 farads 

0 


oscillator solution is 
The following set of initial 


d 2 x 

dt 2 


0 
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CHAPTER HI 

THE LOW-FREQUENCY MAGNETIC FIELDS IN A DISC 
ARISING FROM RADIAL CURRENTS 


Development of the Expression for th e Magnetic Field 
In a metallic disc conductor with free charge, where the permeability, 
M, and conductivity, cr, are isotropic and not time -dependent, Maxwell's 


equations are 



V • D = p 

j 

(8) 

V- B = 0 

5 

(9) 

V x E = -p 

aii 

at 

(10) 


and 


V x B 


M j + M 


d D 

at 


(ii) 


The curl of equation (10) is 


V x (V x E) = - M 


JL 

at 


(V X H) 


( 12 ) 


or 

- - 3 / — 8 E 

V x (V x E) = - M -gj ^ E + € — 

since D = eE, and, by Ohm's law, j = ctE, so that 



15 



Vx (V x E) = - a jx 


9 E 

at 


- e/i 


a 2 e 
at 2 


Now, for a time, periodic field where 


E (q k> t) = E (q fc ) (exp iwt) 


and q^ is a set of orthogonal space coordinates, 


V x [V x E (q k )] (exp iwt) = - i E (q k ) (exp i wt) 


+ € jJLoo 2 E (q k ) (exp iwt) 


(13) 


or 


and 


V x (V x E) 


■ (■-£) 


IjI€<jO Z E 


For aluminum, 

a = 3. 5 x io 7 mho/meter , 


H = (approximately) = 4ir x 10 ' henry /meter , 


«e = (approximately) = — x 10~ 9 farad/meter 

Then, the imaginary term in equation (14) has a coefficient 


e cu 


4. 56 x 10 

w 


18 


( approximately ) 


(14) 


This term is very large compared to unity except for extremely high 
frequencies; therefore, the displacement term, which is real, can be ignored 
at low frequencies. 

Under these conditions, Maxwell's equations become 

V • E = 0 , (15) 

V • i = 0 , (16) 

16 




In an orthogonal coordinate system, the curl of any vector [7] F is 


V x F = 


k (h ’ F ’ J - A 


A_ [4- 

h 3 hj |_9q 3 

,A_ r_9_ 

hih 2 |_9qi 


(hiFj) - 


(h 2 F 2 ) - 


9qi 


9q2 


(h 3 F 3 ) 


(hjFj) 


] 

] • 


The divergence is 
V • F 


i r 

h t h 2 h 3 [c 


(Fjh 2 h 3 ) + (hiF 2 h 3 ) + (hjh 2 F 3 ) . 


8qi 9q2 — — 9q 3 

In a cylindrical coordinate system, the square of the elemental dis- 
tance element is 


ds 2 = dR 2 + R 2 d 9 2 + dz 2 , 


(19) 


and the scale factors are h t = 1, h 2 = R, and h 3 = 1. Now let R = r + r^ so 


that dR = dr, then 


ds 2 = dr 2 + (r+r) 2 d9 2 + dz 2 , (20) 

o 

and hi=i, h 2 = r + r Q , and h 3 = 1. 

The curl then becomes 


V x F 


A 

r 


r + r 

o 


9F 

z 

9(r + r )F 

O C7 

A 

+ 0 

9 F 

r 

9 F 

z 

.99 

9 z 

8 z 

8r 


+ 


A 

Z 


r + r 

o 


9 (r + r ) F 
o 9 

99 


9 F 
r 

9 0 
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and the divergence becomes 


V • F 


r + r 


9(r + r ) F 
o r 

8 r 


r + r 


9F„ 0 F 

— + £ 

0 0 0 Z 


o o 

In the case of a disc with a slowly varying time periodic electric field 
applied in the radial direction, the components of the vectors in Ampere’s law, 
equation (18), are 


. 0 B 9 B 

1 z 0 

r + r 00 ~ 00 ~ M;i r 

o 

0 B 0 B 

- _ £ = o 

0 z 0r ’ 


( 21 ) 


( 22 ) 


and 


i 3(r +y b 9 

r + r 0 r 

o 

The divergence of B is 


1 


r + r 

o 


0 B 
r 

00 


0 


(23) 


r + r 


0(r + r ) B 
o r 

0 r 


r + r 


00 


0 B 
z 

0 Z 


= 0 


(24) 


o o 

(The time periodicity in these equations is not written but is assumed to be 
present. ) 

Since the disc is symmetric about the Z axis, there is no 0 
dependence, and 


B * B( 0) and E * E(0) . 
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From Ohm's law, j = crE for the disc; therefore, if the only nonzero com- 
ponent of j is in the r direction, 

ir ~7 .A 

V x B = r 

From equation (17) , the curl of the electric field is 


_ A 

and since E = E r, 
r 

V x E = V x E r 
r 



3(r E ) 
r_ 

3 z 


3 B 
3 t 


and 

B = B = 0 
r z 

since B has only a 0 component. 
Then 


3 E 3 B 

r_ _ 0 

3 z 3t 

After differentiating, 

3 2 E „ 3B„ 

r_ _ _3_ e_ 

3 z 2 3 1 3 z 

from equation (21). 

From Ohm's law j = 
J r 



crE , so that 
r 


3 2 E 


dz z 


ft 


(25) 


(26) 


20 



The field is applied in the radial direction so that 


E (R, z, t) = E (R, z) (exp i cot); (27) 

r r 

then 

8 2 E 

5 — = ionru E (R, z) (exp i cot) 

3z^ r 

From equation (15), 

V- E = 0 , 


or 

8(RE ) 

R — " ° ' 

and, from equation (27) , 

[RapE (R, z) (exp i cot)] = 0 
K d Jr! r 


This implies that 

RE (R,z) = g(z) (exp i cot) . (28) 

r 

The function g(z) must be an even function since the electric field 
must have the same magnitude on both sides of the disc. 
g(z) = g(-z) 

The maximum electric field on the disc occurs when exp i co t = 1 
and when the variables R and z have the values 
R = r 

o 

z = ±z 


max 

The function g(z) 


g(z ) = r E 
& max o or 


will then have the value 


Then, from equation (28), 


E r (R,z, t) = 


(exp i wt) 


It has been shown that 


r 8 , _ . 

T~ = TT (<rnE) ; 


‘ R"k [/^dT 1 ] ^P 1 " 11 = ~t <«piwt) 


(^E r ) = - lC ^ afX g(z) (exp i co t) 


Then 


d 2 g(z) 

dz 2 


= i co crp g (z ) 


(D 2 - i co op) g(z) = 0 


This equation has a solution of the form 

g(z) = Cj(exp sTiio ajl z) + C 2 ^exp(- 
Since the function g(z) is an even function, 


g(z) = g(-z) 


I i co a p z 


Cj - C 2 , 


then g(z) will have the form 


g(z) = C 1 j (expN/ico up z 


) + £ exp(w ico op) z | 
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At 


so that 


or 


Since 


then 


z = ±z 


max 


g(z ) = r E 
max o or 


C, = 


r E 
o or 


exp 


^ i oo a n z + exp (- \] i oo crfi ) z 
max x 7 


max 


g(z) = r Q ^ or j(exp'\/ia) cr fi z) + £exp (- \lico ufi ) zj | 


-r | ^expN/ia) crju z max ^ + exp^- n/ i u> cr/u ^ 


max 


l = exp 


1 7 r 
2 ’ 


rr- 1 7T 7T . . 7T N/1T , . 

vi = exp — — = cos — + i sin — = — — (1 + i) 

4 4 4 2 


(31) 


/ v ( \l2afJuo \ ( W2 crju^ \ 

s<z) * r o E or l exp ~~2 z ) V eXP “ 2 z / 

/ W 2 cr/u an 

r - * — > 


\IYajjuo \ 

exp ^ j z 


exp - 


* ((' 


'J2<rn ^ 

exp_ 5 z max ) V exp 2 


1^2(7^00 


exp 


iax ) 


max 



After expanding, and letting the second term be A , 


g(z) = r E 
6 o or 


(« 


n/ 2 crp co . n/ 2 cr/u co \ f \/2afi u 

i cos z - sin z M cos z 


. . , , to \ /. 

+ 1 sin — z J + l 1 cos 


. n/ 2 op < 

2 


n/2 ct/xco . \/2 op 

z + sin — 

2 2 


-0 


( 


x \ cos 


n/ 2 ap o) 


.« ) 

z - i sin z 1 


. n/ 2 up < 

2 


-r A 


or collecting terms and writing out A , 


= r E (c 
o or \ 


n/ 2 <7jU CO 


2 ^/2 (T)U CO 


g(z) = r E \ cos" — z - sim r" z 


*( 


2 V2 (S\X CO 


. o v2otuo 

cos" — — "T z - sin z 

2 max 2 max 


) 


z is the largest z dimension of the disc and is equal to 
max 

3. 175 x 10 -6 meters. For aluminum, a{x = 44 (approximately); then, for a 
frequency of 13 172 hertz, ^2 co _ 5. 38 x 10 2 . 


\l2 (JLi CO ^ , \To-Jujo~ 

Under these conditions, sin z = 0 and cos — z = 1 

z z 

approximately. Therefore, g(z) is constant over the dimensions of the 
disc and the skin effect may be neglected at this frequency. 

It has been shown that 


E (R, z, t) = 
r 


g(z) 


R 


exp i co t 


and 


B = 0 


24 
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so that 


R 9E 0 


from equation (24). 

Then RB 0 =f(-z,t) +C, utmost. 
From equation (25) , 

9 B 


at 


r at + ~ r at f(z,t) 


Then, 


a E 


= J_ d g(z) 


(expicot) = - ~ ~ f(z,t) , 


a z r dz * w " "R3t 

/ a"t dt = - — ^ z “~ j (exp icot) dt , 

f(z,t) = ~ (expicut) + C 


B„ = 


0 coR 


dg(z) 
d z 


(exp ico t) + C 


and 


B e - 7 ^r r 1 [expi(cot+ f + C r )] 


From equation (21) and Ohm's law, 
aB e 


Sz = - af,E r 


and 


^6 
d z 


1_ d 2 g(z) 


co R 


d z 1 


jexp i (cot + ~ +c )] 


• (32) 


(33) 
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Now, 


and 


/, „ ( -> ^2(7f«tD . , \/2 (TJU co \ 

g(z) = r E loos' 1 — z - snr z J 

o or \ 2 2 / 


♦( 


2 V2 CrfJLGO . 2 ^2 ( THU) 

cos 6 z - snr z 

2 max 2 max 


) . 


dg(z) „ ( rz . n/TotTw \JlTajjuo 1 

— “ = - r E \ \l2cr nco sin z cos z } 

dz o or \ 2 2 / 


*( 


„ 2 n/2 c<; 


„ . , \/2crpw 

cos‘ z - snr — — z 

2 max 2 max 


d 2 g(z) „ ( 2. ^2 crfiw . 9 \T2an 

— f - o = r E <jh<jo[ cos* — - £ - — z - snr — 

dz 1 o or ^ \ 2 


E 0 


*( 


o \p2au U) . 9 n/ 2 cru cu 

cos^ — c — z - snr z 

2 max 2 max 


d 2 g(z) 

dz 2 


= q/uw g(z ) 


(34) 


and by substituting for E 


c ru 

- qjuE^ = - g(z) exp iwt 


Then 


7 T T 

an g(z) exp iw t = an g(z) exp i(wt + - + C ) 


and for this condition to be true 


f 
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Then, from equation (33) , 
E 


B 


— r 

( 


ylz (j u, co / 

or ( 


to R 


^2 au co \J 2 cru co 

sin — ^ — z cos — - c - — z 

z z 


) 


♦ ( 


cos 6 


\/g (JJX 0) 


max 


- sin' 


2 n /2 CTfJL co 


max 


3 


exp i co t 


(35) 


Therefore, the magnetic fields in the disc are in time phase with the 

electric fields, and, since the cosine is an even function and the sine is an odd 

function, the direction of the magnetic field reverses between z = z 

’ ° max 

and z = -z 

max 


The Magnetic and Thermal Pressures on the Disc 

The direction of the magnetic fields in the disc are such that if the 
current flows from the center of the disc to the rim, the field on the disc 
surface at the muzzle end of the plasma generator will be directed clockwise 
when viewed from the muzzle. Conversely, the field on the opposing surface 
will be counterclockwise. Thus, the force on the plasma will be directed into 
the disc on both sides of the disc. 

The magnetic force will be: 

j x B = - — — B x (V x B) , 

M 

j x B = - —V B • B + — (B • V) B 

fxB = - — + — B • V B . (36) 

J ix 2 y. 
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g2 

The term is called the magnetic pressure, and the second term can be 

considered as force directed parallel to the magnetic field [ 8] . 

If the gas dynamic pressure exceeds the magnetic pressure, the disc 
shape of the plasma may be distorted. It has been shown [equation (35)] that 


B 


T o E r>r ^ . ^apco n/2j 


r 

- 


® |_ co R 


(T IX CO 

sin z cos z 

2 2 


) 


*( 


o \/ 2 auco . 2 \/ 2 au co 

cos'* — z - sin — — ^ — z 

2 max 2 max 


\] 

f exp i co t 


and, if, 


n/2o~p 


co 


z << 1 , 


then 


cos 
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„ _ 


z = 


„ , . \ 1 2 a/i CO 

i and sin r ^ — z 


V2 


crp co 


z (approximately) 


Under these conditions, 


au r E 

o or . 

•r = — — z exp i co t 

e R 


(37) 


Since the resistance of the coaxial gun is negligible, it may be considered 
as a coaxial capacitor made up of the two conductors which are then short- 
circuited by the aluminum disc. Then, the inner conductor may be considered 
a line charge. The field outside the radius r^ is given by Gauss' law for a 
cylinder where 

f e E • dS = Q 
J o 
s 
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and 


e E (2irr )i = Q , 
o or o 


or 


q t 


E 


or 2 7r e r 
o o 


where Q is the charge density per unit length. 

Using this value for E will yield the following equation for magnetic 
field strength: 



Now , 

E = -V V 

and, for a cylindrical conductor of line charge density Q L> 
Qr 


E = 


r 2ire R 
o 


dV 

dr 


or 


V_ - V 
R ro 

max 


Q t R 

L , max 

In 


2 7T € 


and 


Q L ” 


2ir e V 
o ro 


R 


In 


max 
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when V 


R 


0 . Then, the magnetic field is 


max 


o\x V 


B * = 


ro 


R 


z exp i co t 


R In 


max 


The maximum magnetic field strength will occur when 


z = z 


max 


R = r 

o 

and 

exp i co t = 1 

Then, with the voltages and dimensions used in the experiment, the maximum 
magnetic field will be, assuming the full voltage is applied to the disc, 

(max) = 3.933 x 10 2 weber/ square meter at R = 

and 

^B(min) = 65.73 weber/ square meter at R = R 
0 max 

The magnetic pressure term is 

10 

p = 6. 16 x 10 Newtons/ square meter at R = r 
m o 

and 

9 

p = 1.72 x 10 Newtons/ square meter at R = R 

max 
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If the neutral gas and plasma were in thermal equilibrium, the disc shape 
of the plasma might be distorted. There is no evidence of such distortion, 
and to the contrary the plume leaving the muzzle is pointed as it would be 
under the influences of the magnetic field arising from current in the center 
electrode. This behavior is depicted in Fig. 6. The pointed plume is 
caused by the difference in magnetic pressure near the center electrode 
and the magnetic pressure near the outer electrode. 

The ionized aluminum cannot expand into the accelerator since it is 
under the influence of the j x B forces arising from the magnetic field 
around the center electrode and the current in the plasma. However, the 
gaseous aluminum will expand into the accelerator since it is not acted upon 
by the magnetic field. Also, because of equipment limitations in available 
vacuum pump speeds, the atmospheric leakage into the range made it impos- 
sible to achieve a vacuum below 0. i Newton/ square meter. This much 
atmosphere within the coaxial gun provides an easy electrical breakdown path 
between the electrodes. The net result is that instead of the idealized physical 
situation considered in the development of the theory a much more complex 
situation actually exists within the plasma gun. After a very short time the 
gun is filled with electrically conducting gas and acts only as a resistive 
element in the circuit, except for that portion of the gun which is occupied by 
the insulating plug ( Fig. 7 ) . 


3i 
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FIGURE 6. PLASMA PLUME LEAVING THE MUZZLE 
OF THE COAXIAL ACCELERATOR 






CHAPTER IV 


DISCREPANCIES BETWEEN THEORETICAL CALCULATIONS 
AND EXPERIMENTAL RESULTS 

By examining Figures 8 and 9 it can be seen that a large frequency 
discrepancy exists between the calculated current versus time curve and the 
experimental curve. Also, the difference between the calculated velocities 
and the velocities obtained experimentally is significant. The simple theory 
used in the calculations is evidently not sufficient to account for the behavior 
of the system. Discrepancies have also been noted in the literature [3]. An 
explanation is offered which may account for the differences between the 
theory and the experiment. 

The aluminum disc used to load the gun weighs 39 milligrams, while 
the calculated mass of material to be ionized is 2. 003 milligrarfis when a 
value k = 1 is used in the expression 

V 2 C 2 L 2 
o 

m ~ 2k(L + Ax ) 
o o 

When the aluminum disc is vaporized, the excess mass expands into 
the coaxial gun and, because of collision processes between aluminum atoms, 
provides ions to spread the discharge in a direction opposite to the accelerating 
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FIGURE 8. THEORETICAL CURRENT VERSUS TIME CURVE 



FIGURE 9. CURRENT VERSUS TIME WHEN FIRING COAXIAL GUN 
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j x B forces. To augment this spreading discharge, a considerable amount 
of atmosphere remains in the gun at a pressure of 0. 1 Newton/square meter. 
This gas is sufficient to cause a breakdown between the gun electrodes when 
a large difference in potential is applied. 

That portion of the coaxial gun which is involved in the discharge now 
acts as a resistive element in the circuit. Current must still flow in the two 
electrodes but the current per unit length of electrode is no longer constant. 
Thus, the magnetic field produced by the current in the inner electrode is 
changed and the inductance of the circuit approaches the value L^. 


36 



I 


CHAPTER V 

THE ENERGY AND MOMENTUM OF THE PLASMA 

The Energy Contained in the Plasma 
If the previous assumptions are correct, the experimental discharge 
current curve can be explained. The time interval between the initial current 
value and the first zero in current is longer than the time interval between 
succeeding zeroes of current. The time interval between successive zeroes, 
after the first, agrees with the period of the ringing circuit without the gun. 

The mass of the ions accelerated by the magnetic field cannot be 
obtained from the simple theoretical solution, but must be approximated from 
one of the three experimental curves for current and velocity. They are the 
ringing frequency current Versus time curve for the circuit without the 
coaxial gun (Fig. 2) , the current versus time curve for the circuit when 
firing the coaxial gun (Fig. 9) , and the distance versus time curve for the 
advancing plasma plume in the vacuum chamber (Fig. 10) . 

The equations for current, frequency, and logarithmic decrement 
shown in Appendix B are 
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The current expression when solved for the current at the first maxi- 
mum produces an expression which is proportional to the original voltage 
stored in the capacitor. Without changing any circuit elements, but with a 

voltage applied in the opposite sense and of a value V* , the second extrenum 

o 

in the discharge current curve could be obtained after the same time interval 
38 



used in obtaining i . . The ratio of the absolute values of the first two 

1 max 

current extremums is then equal to the ratio of voltages on the capacitor 
giving rise to the current extremes. 


1 max 


2 max 


V ! 

o 


Then, 


v' = . . 

o I 1 


max 


i max 

Since the energy in the capacitor bank is E = C V 2 , the energies are 

A I 

E 1 = -j C V q 2 = 34. 56 kilojoules and E 2 = ~ C W 2 = 14. 8 kilojoules. The 
energy difference between these two is the energy dissipated in the circuit 
during the first half-cycle of the current discharge curve. 

The distributed resistance of the circuit without the coaxial gun has 
been calculated from the ringing frequency current discharge curve , and 
the resistance of the coaxial gun was measured using a Wheatstone bridge. 
The distributed resistance of the circuit, including the coaxial gun, can be 
calculated from the discharge current curve and is 


= 1. 87 x 10 -2 ohms 


The resistance of the plasma is 

R p = R t - R Q - Ri(0.67) 

Rp = 8. 819 x 10 -3 ohms 
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Then, the energy dissipated is the plasma is 


E p = AE 



where 


AE = Ej - E 2 and E p = 9. 3 kilojoules 

By examining Figure 10 it can be seen that the plume front leaves the 
muzzle of the coaxial gun after 12. 5 x 10 -6 seconds. A Simpson's rule of 
integration of the data from the experimental current curve where 


t 

E = / I ! R dt over the first half-cycle of the discharge and over the first 
o 

12. 5 x 10 -6 seconds of discharge can be used to calculate the fraction of the 
energy absorbed during the first 12. 5 microseconds of the first half-cycle. 


f 


12. 5 x 10“ 6 

/ 


o 


34 x lO" 6 

/ 


o 


I 2 R p dt 


I 2 R p dt 


0. 32 


Then, 

A Ep = 0. 32 x 9. 3 x 10 3 = 2. 98 x 10 3 joules . 

The fraction of this energy which is actually retained by the aluminum 
gas is unknown because the heat exchanged between the gun and the electrodes 
has not been measured or accounted for. There are certain observed 
phenomena that indicate that a large amount of heat is delivered to the 
electrodes. For example, after the discharge, the center electrode is 
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slightly roughened, as though material had been removed, and a noticeable 
amount of shiny metal, assumed to be aluminum, has been deposited over 
this. The outer cylindrical electrode is coated on the surface with a grayish 
black material, assumed to be aluminum oxide. There is no quantitative way 
of measuring the heat absorbed from the gas. Also, no attempt has been 
made to distribute the energy among the terms in the energy-balance equation 
for the partially ionized aluminum gas. 

As a first approximation, consider all of the energy absorbed by the 
plasma as being converted to kinetic energy of the plasma. Then, 


The quantity v , the average velocity of the plasma, is obtained from 

clV 

Figure 10 which shows the time needed for the plasma plume to travel from 
the muzzle of the gun to a point 25. 4 centimeters downstream. 
v &v = 23 kilometers/second 

Then, 

2 x 2. 98 x IQ 3 
m “ (23) 2 x 10 6 

m = 1. 13 x 10 -5 kilogram , 

m = 1 1. 3 milligram 

Since only the ions are accelerated by the magnetic field, and if it is assumed 
that all of the aluminum is vaporized, the percent ionization is 
11 . 3 

r] = — — — = 29 percent 
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The ions are being accelerated and it is assumed that they fill the gun 
from the original position of the foil out to the muzzle. 

Unquestionably, these values for the energy and mass are too high 
since the energies of ionization, excitation, particle interaction, degeneracy 
(Fermi correction) , radiation, and plasma oscillation [9, 10, 11] are not 
taken into account and the heat lost to the electrodes is ignored. 

Transfer of the Plasma Momentum to a Projectile 
Since the measured velocity is an average velocity, it will be used 
to calculate the momentum transfer from the ionized gas to the simulated 
meteoroid material. It is assumed that the number of ions contained in any 
cylindrical shell concentric about the center electrode and of radius less than 
the radius of the outer electrode is constant. This must be true since current 
is flowing radially between the two electrodes. 

Now consider a glass bead 5 x 10 -5 meters in diameter to be placed a 
distance r from the axis of the gun and at the muzzle. Since the difference 
between the o. d. of the center electrode and the i. d. of the outer electrode is 
3. 235 centimeters, a total mass of 1. 75 x 10“ 8 kilograms resides in the ring. 
The mass intercepted by the glass bead will be 

Am = —A L 

27rr x 5 x 10“ 5 

At the surface of the center electrode, this mass will be 
Am = 1. 72 x 10“ 11 kilograms 
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This is true assuming a laminar flow so that the geometric cross section of 
the glass bead may be used in calculating the intercepted mass. 

Now , if the collision process between the aluminum ions and the glass 
bead is elastic so that both linear momentum and kinetic energy are con- 
served, the velocity of the glass bead may be calculated. 

The glass bead has a density of 2. 51 and a mass of 1. 64 x iO -9 
kilograms: 

(Am) vj = (Am) V 2 + m^ V 

-|-(Am) vi 2 = y (Am) v 2 2 + -j- m Q V 2 

This set of equations must be solved for V, the velocity of the glass bead, 
when vj = 23 kilometers/second: 

V = 0. 48 kilometer/second 

No provision has been made for the momentum transfer from drag forces 
between the plasma and the projectile. It has been considered negligible in 
performing this calculation. 

Also the collision process between plasma and projectile is not com- 
pletely elastic since the projectile is ablated to some degree. The inelastic 
collision will result in a slower velocity projectile, however, so that the 
elastic collision calculation is the most favorable from the point of view of 
determining maximum possible particle velocity. 

It has been determined that the coaxial plasma accelerator is not 
a suitable device to propel projectiles to velocities of the order of 
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20 kilometers/second unless the accelerator is modified in a way which will 
increase the plasma density. 

Several techniques have been tried for increasing the plasma density 
(e. g. , converging nozzlesof dielectric and conducting materials). 

When a stainless steel nozzle was used, no beads striking the target 
were detected. It was noted that spherules of condensed aluminum, 1 to 2 
millimeters in diameter, appeared on the wall of the nozzle after the gun 
was fired. Since it appeared that too much energy was absorbed from the 
plasma by the stainless steel nozzle walls, a dielectric material was tried. 

A nozzle was constructed with a commercially available material 
called "Grade-A Lava. " When this nozzle was used, velocities of 
2. 3 kilometers/second were attained with 1. 05 x 10 -2 -centimeters diameter 
glass beads. Although the walls of the nozzle were blackened, no noticeable 
deposit of aluminum on the nozzle wall was detected. However, since the 
velocity attainable was only about 10 percent of the desired velocity, it was 
necessary to find a different way to increase the plasma density. 

A coil wound in a conical shape, concentric about the longitudinal axis 
of the gun and energized by the plasma itself, has accelerated glass beads of 
0. 66-millimeter diameter to velocities of 7 to 8 kilometers/second and beads 
of smaller diameter to still higher velocities. The theory of the coil has not 
yet been fully developed but is under investigation. 
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CHAPTER VI 


SUMMARY AND CONCLUSIONS 
Summary 

It has been shown that the simple circuit theory of the coaxial plasma 
generator is not sufficient to adequately explain the physical phenomena which 
occur when the plasma generator is loaded with a disc of conducting material 
and then energized. In particular it has been shown that the magnetic fields 
developed by the current flowing in the disc are not of sufficient strength to 
keep the disc from exploding. 

It has also been shown that this coaxial plasma generator is not suitable 
for accelerating projectiles of a useful size range to velocities of the order 
of 20 kilometers/second unless some means of increasing the plasma density 
is used. 

It has also been experimentally determined that when the density of the 
plasma is increased, projectiles of a useful size range and velocity have been 
produced. 
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Conclusions 


By utilizing the coaxial plasma generator and a self-energized plasma- 
containment coil, a useful device may be designed for accelerating small 
projectiles to hypervelocities in the range of 20 kilometers/second. 

The theoretical understanding of this system and the parameters to be 
varied to obtain predicted results have yet to be discovered. 



APPENDIX A 


INDUCTANCE AND CAPACITANCE OF COAXIAL CONDUCTORS 

For a hollow conducting cylinder carrying a current the magnetic field 
at the points P and P' of Figure A-l can be calculated from the Biot- 
Savartlaw, where 



FIGURE A-l. HOLLOW CONDUCTING CYLINDER 



Referring to Figure A-l, points P and P* lie on a diameter and are 
equally distant from the longitudinal axis of the cylinder. Consider pairs of 
current elements, Idij and Idi 2 > aQ d Idij and Id ^2 where I is the 
current per unit length of perimeter of the cylinder, and use the Biot-Savart 
equation to compute the resulting fields 


dB = 



df 2 X W 2 
w 2 3 


and 


dB 


4 7 r 


dii x w! 


W| 


»3 


d l\ x- 1 


w 2 


W 2 


' 3 


Since the magnitudes of the df and w terms are such that 
df ! = df 2 = df / = di 2 ' an d Wj = w 2 = w/ = w 2 ' , 

the only nonzero components of dB and dB' are dB^ and dB^ , and since 
these are components of oppositely directed vectors along the diameter, the 
vector sum is 


dB = r dB - r dBJ = 0 . 

K x\ 

Therefore, if all current elements are calculated for every pair of 
interior points of the cylinder, the resultant magnetic field will be zero. 
There is no magnetic field within a hollow-conducting cylinder carrying a 
longitudinal current. 
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The field surrounding the center conductor of a coaxial line may be 
found using the circuital form of Ampere's law, where 

<f B • dT= n I 
J o 

for a path linking the current. For a path tangential to the field 
f B • di = B f di 

and 

B di = 2ir rB 
for a circular path. 

Then, 

ju I = 27T B r , 
o 


or 


B = 


V: 

27rr 


These results may now be used to find the inductance of a coaxial pair 
of conductors. The total magnetic flux linkage for a length Z is 
b y. I b 


A = z / B dr = z / 

a a 


dr 


z u 

, o , b 

A = — In — 

2 7r a 


The inductance L is defined as 




M z . 

o , b 

— In — 

2n a 


and the inductance per unit length is 
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To calculate the capacitance for the coaxial conductors , start with 
Maxwell’s equation, where 
V • D = p 


V • E = 


= 


From the integral form of Gauss' law, 


J E • ds = f — d t = 


Q t 


where the charge Q is a line charge 


Qi 


E r (27rr) = — 


density, or 


since 


J ds = 27rr 

per unit length. Then, 

Q 

„ L 

E = 

r 2 7r € r 

o 


Since 


E = - VV 


and 


E 

r 


dV 

dr 
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or 


R 


- 1 


max 


R 


max 


dV = / E dr , 


and, if 


= 0 


max 


then 


R 


/ 


max Q 


o r 


2tt e 


dr 

r 


Q t R 

L , max 
In 


2ir e 


By definition, 


Q 

c = ■ 

C V 


so that the capacitance of the coaxial cable is 
2ff e 


C = 


R 


farad/meter . 


In 


max 
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APPENDIX B 


LCR CIRCUIT THEORY 


For the circuit consisting of an inductance L, a capacitance C, and 
a resistance R connected in series, the differential equation describing the 
circuit is 


L f • * 1R+ ■ 


dq 

Now i = — , so that 
at 


t d2 q , R i^ + i 

dt 2 dt C 


= 0 , 


and, by differentiation, 


r d2 i D di i 
L dt 5 ” + R dt + C 


= 0 


then 


t , b ^a, 1 

L ^t 7 dt C 


= 0 


has a solution of the form Q = A (exp yt) . Then, substituting into the dif- 
ferential equation, 


L y l + Ry + — = 0 , 
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so that 


= _ _R / R 2 1 

“ 2L ± n/ 4l7 " LC 


Letting 


a s 


R 

2L 


and 


p 


s ± /-5| - -L 

n/ 41? LC 


then 


7= -« ± 0 , 


so 


q = (exp - at) [A a (exp /St) + A 2 (exp - /St)] 

and 


i = (exp - at) [ (/3 - a) A t (exp /St) - (/3 + a) A 2 (exp - /St)] 
Now q = q Q and i = 0 at t = 0, or Q () = A 1 + A 2 , and 
(/3 - QjJAj = (/3+ a) A 2 , 


A 


a -/3 
2/3 


and 


A - a 
Al " 2/3 
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Then, 


q= — (exp - at) [(a + /3) (exp fit) -(a -fi) (exp -/3t)] 


If /3 is complex, then 


R 2 1 

"41? " LC < ° ’ 


or 


r!<4 £ ■ 

and a periodic solution is possible. 

If a solution is assumed for q» where 

q = A q (exp - at) sin (w^ + <p) , 

and is substituted into the differential equation 

[ ( a 2 -a> 2 ) L - a R + ^ sin( c 0 t + <p) + ( -2 a co L +a> R) cos( u t+cp) = 0 . 
O LI O O O O 


Then 



Using the initial conditions 


q = A sin cp 
o o r 


and 


cos cp - a sin cp - 0 , 

then 


tan cp - 


c o 

o 

a 


A = 
o 



sin cp 
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or 


CO 2 + a 

i = y ° 


q (exp - ck ti . [ oj cos <07 t + cp) -a sin(o> t + cp)\ 
0)0 o o o 

o 


By expanding, 


l = - 


o) 2 + a 2 

w ' q Q (exp - at) sin o; q t , 


or 


i = i (exp - at) sin u> t , 
o K o 


where 


i = 
o 


co 2 + a 2 


a; / o 
o 


and 


= ^o _ _L_ I 1 _ 

'’o 27T 27T v LC 


R 2 

iU 


To find the time when the current has an extremum value, differen- 
tiate the expression 

j 2 
o 


' w 2 + a 2 


i = - 


co 


q (exp - at) sin co t 
o o 


is differentiated and the derivative set equal to zero. Then, 

co 

. , o 

tan co t = — 
o max a 

Since there are recurring maxima and minima, the time interval between 
successive maxima is 


co (t + T ) = m t +27 r , 

o max o o max 
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where 



If i* and i M are two successive maxima of current, then 


i_^ _ (exp - at) 
i" (exp[-Q!t - arT^] 


exp a T 

o 


and 

In 
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APPENDIX C 


COMPUTER PROGRAM 


A fourth -order Runge-Kutta integration scheme written in MARVES 
language [12] was used to solve the coupled circuit equations 


< L „ + “> 5F + ( R o +R < x+R ^ x f + 5 = ° 


dx \ dV 


and 

d 2 x L o + ^o / dV\ 2 

dt 1 V \ dt / 

o 


The variable and constant names used in the program are: 


L = LO 
o 

A = LAMBDA * 


x = X 



= DDVT 


R = RO 
o 

R! = R1 
R 2 = R2 

dx 

dF =DXT 
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d 2 x 


~w 


= DDXT 


dV 

dt 


= DVT 


V = V 

c = c 

t = TIME 


The computed current was 


CUR = - (C) (DVT) . 


The flow chart for the main program is shown in Figure C-i. The 
subroutine SETONE flow chart is shown in Figure C-2. Within the SETONE 
subroutine is a Simpson's rule integration of the quantity I 2 R 2 which provides 
a calculated energy dissipation for the plasma. The integration is carried 
.out over a time period equal to the duration of the first half -cycle of the cur- 
rent discharge curve ( Fig. 9) . 

The Simpson's rule integration is called SIMPSO and is called out in 
subroutine SETONE as SIMPSO. The SIMPSO subprogram flow chart is shown 
in Figure C-3. The flow chart for subroutine SETT WO is shown in Figure 
C-4. The computer program listing is also included. 
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MAIN PROGRAM 

COMMENT: LANGUAGE FORTRAN 



FIGURE C-l. COMPUTER MAIN PROGRAM 
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SUBROUTINE SETONE 

COMMENT: DIFFERENTIAL EQUATIONS FOR X *0.67 AND 
TIME -200 At 
LANGUAGE MARVES 



FIGURE C-2. COMPUTER SETONE PROGRAM 
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COMMENT: SIMPSON'S RULE 

LANGUAGE FORTRAN 



FIGURE C-3. COMPUTER SIMPSO PROGRAM 
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SUBROUTINE SETTWO 

COMMENT: DIFFERENTIAL EQUATIONS FOR X > 0.67 AND 

TIME -200 M sec 
LANGUAGE MARVES I 



FIGURE C-4. COMPUTER SETTWO PROGRAM 
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00101 

1 * 


SEAL LO.L AMBDA 


DD1D3 

2* 


COMMON/ BL0CK1/L 0 . 1 A MBD A » RO » R 1 » R 2 »C . C ONS T/B l OCK 2 /V »D V T»D D VT 

t X*0 X T • 

00103 

3* 


SODXT .TlME.DTST EP 


00 ID 4 

4* 


DIMENSION CUR1 ( 400) .TImEI (400) .FlDX< 12) .FLO Y( 12) tFLDZf 12) 


00105 

5* 


COMMON/ BL OCK 3/ CUR 1 *TIM El. NUMBER 


00106 

6* 


CALL I CENT ( 1 ) 


00107 

7* 


READ( 5. 1) Fl OX. F. DY. FLD2 


00125 

8* 

1 

FORMAT (12A6/12A6/12A6) 


00126 

9* 


cau setone 


00127 

ID* 


CALL SETTWO 


00130 

11 * 


00 10 I : 1. NUMBER 


00133 

12* 

10 

CURKI) z - C U R 1 ( I ) 


00135 

13* 


CALL QJ IK 3L (- 1 . U . . 200 .• - 3 . . 3 . . 1 H* . FL DX . FL DY NU MBER. T I«E1 . 

CUR1 ) 

00136 

14* 


DO 2 I : 1. NUMBER 


00141 

15* 

2 

CURKI) = CUR 1 ( I ) * 1 • E+5 



00143 

16* 

CALL 

GUIK3L (-1.0..200, . -3.E5 .3.E5. 1H*.FLDX.FL0Z.-NUMBER*TIME1.CUR 

00143 

17* 

S 1 ) 


00144 

1 8* 

CAL L 

ENDJOB 

00145 

19* 

STOP 


00146 

20* 

END 



END o r compilation: no diagnostics. 
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CO 



O 


1 * 
2 * 

3 * 

4 * 
5 * 
6 * 
7* 
8 * 
9* 

10 * 
11 * 
12 * 
13* 
l 4 * 
15* 
1 6 * 
17* 
18* 
19* 
20 * 
21 * 
22 * 
23* 
24* 
25* 
25* 
27* 
28* 
29* 
30* 
31 • 
32* 
33* 
34* 
35* 
38* 
37* 
38* 
39* 
40* 
41* 
42* 
43* 
44* 
45 * 
46 * 
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SUBROUTINE SETONE 
REAL L 0* l AMBDA 

COMMON/BLOCKl/t 0 » l A MBOA * RO * R 1 • R2 * C * CON ST /B t OCK 2/ V *0 VT tDOVT.X *OXT# 
SDOXTf TIMEtOTSTEP 

0 IHENSION CUR1( 400) #FCUR< 200 ) .TImEI (400) 

CQMMQN/Bt OCX 3/ OJR1 • TIM El .NUMBER 

INITIALIZE 

NUMBER* ~ Q 

INTEGRATION ACCURACY r 16 
SKIP PAGE 
3 CONTINUE 

I NPUI NAMEL 1ST (5*N AM El )l 0»L AH8DAt7 Dt R1 1 R2* C» V »X # DVT » OXT • TIME* DTSTE 
IP 

OT STE° = 0TSTEP/10. 

6 CONTINUE 

OUTPUT TITLE RUNGE KUTTA INTEGRATION 
SKIP LINE 

CONST - (LO ♦ L AMBDA*X ) / ( V**2*l AM80A) 

OUTPUT TITLE DEFINED AND CALCULATED CONSTANTS 
SKIP L INE 

OUTPUT LIST L0*LAMBDAtR0tRlfR2/C*C0NST 
SKIP L INE 

OUTPUT TITLE INITIAL CONDITIONS 
SKIP LINE 

OUTPUT LIST TI ME *DXT * X tD VT • V 
SKIP L INE 

differential equations 

2 DOXT r CONST *DVT **2 

OOVT =: - ( ( R 2 ♦ RU ♦ R1*X ♦ t Av&DA *D X T ) *D VT ♦ V/C)/(lU ♦ l AMBDA *X ) 

1 CONTINUE 

integration 

METHOD RUNGE KUTTA ORDER - 2 * X t OX T t ODX T ♦ V ♦ OVT t DDVT 
EVENTS 
Cl = OXT 
CUR = -C*DVT 

OUTPUT LIST TlME.VtCUR»DXT»X 
NUMBER - NUMBER * 1 
FCUR( NUMBER) = CUR**2*R2 
CURl(NUMBER) = CUR/1. E+5 
TI ME1( NUM BER) □ TImE* 1.E*6 
I F ( X *GE • 0 . 6 7 ) GO TO 7 
EV ENT f EX IT) T IM E = 2.E-3 
TERM INAL COMPUTATIONS 

7 AREA z SlMPS0(OTSTEP*NUMBEtffFCURJ 
SKIP PAGE 

OUTPUT LIST AREA 

return 

ENO 


OOlOl 

1 * 


FUNCTION SIMPSOl HB.N .Y ) 

0010 3 

2* 


DIMENSION Y(N) 

00104 

3* 


N2 - N/2 

00105 

4* 


N3 : ( N ♦ 1 )/2 

00105 

5* 


IF( N2.NE.N3) GO TO 19 

00110 

6* 


AREA z (Y(N> - Y ( N - l)»*HB/2 

00111 

7* 


N - N - 1 

00112 

8* 


60 TO 20 

00113 

9* 

19 

AREA = 0. 

00114 

10* 

20 

XI T Y ( 1) ♦ Y ( N ) 

00115 

11* 


X 2 : 0. 

00116 

12* 


X 3 r 0. 

00117 

13* 


Nl - N - 1 

00120 

14* 


DO 1 I : 2 » Nl » 2 

00123 

15* 

1 

X 2 = X2 * Y ( I) 

00125 

16* 


N2 = Nl - 1 

00126 

17* 


DO 2 I : 3.N2.2 

00131 

18* 

2 

X3 z X 3 * Y 1 1 » 

00133 

19* 


SIMPSO = (HB/3.1 *<X1 * 4.*X2 

00134 

20* 


return 

00135 

21* 


END 


AREA 


END OF COMPIl ATION: 


NO diagnostics. 


Ci 

oi 


05 

05 


1* SUBROUTINE SETTWO 

2* REAL LOtlAHBDA 

3* COMMON/ BL OCX 1/L 0»L AM BDAfRO»Rl »R2»C»C0NST/8L 0CK2/V* DVT »ODV T«X. OX T» 

4* JDDXT »TI ME iD TSTEP 

5* DIMENSION C'J R 1 < 400 ) t T IM El < 400 ) 

6* C 0MM0N/B10CK3/CUR1»TIME1 »NU«BER 

7* initialize 

8* DTSTEP = DTSTEP*1U. 

9* INTEGRATION ACCURACY = 16 

ID* SKIP PAGE 

11* OUTPUT LIST LO.L AMBDAiR0tRl.R2/C. CONST 

12* SKIP LINE 

13* OUTPUT LIST T IM E # OV T » V • OX T . X 

14* SKIP LINE 

15* DIFFERENTIAL EQUATIONS 

16* D DVT = -MRU ♦ R2 + R1*0.67)*DVT ♦ V VC > / ( L 0 ♦ 0.67*LAMBDA) 

17* INTEGRATION 

18* METHOD RUN6E KUTTA ORDER - 2 • V *0 VT tDDV T 

19* EVENTS 

20* CUR r -C*DVT 

21* OUTPUT LIST TIME.V»CUR 

22* NUMBER = NUMBER * 1 

23* TI ME1( NUM 8ER) ^TlME*l.£ + 6 

24* CUR1 (NUMBER )= CUR/1. E+5 

25* EVENK EXimiME = 2.E-4 

26* TERMINAL COMPUTATIONS 

27* RETURN 

28* END 



APPENDIX D 


BEHAVIOR OF AN IONIZED DISC CONDUCTOR IN A COAXIAL 
PLASMA ACCELERATOR 

The magnetic fields within the ionized disc have been calculated and 
shown to be clockwise when viewed facing the disc. The fields are directed 
oppositely on the two sides of the disc. 

The magnetic pressure term arising from these magnetic fields has 
been calculated for the particular conditions described in the body of this 
report. The pressure has been shown to be 36 times as high near the center 
electrode as it is near the outer electrode. The pressure is caused by the 
magnetic fields within the disc which create a force inward on the current 
on both sides of the disc. 

The magnetic field created by the current flowing in the center 
electrode of the coaxial accelerator toward the disc creates a magnetic field 
around the center electrode in the same direction as the field on the inner 
disc face caused by the current in the disc. The field of the center electrode 
has a magnitude of (Appendix A) 

Ml 

T> = ° 

27 rr 
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The charged particles are acted upon by the j x B forces arising from the 
current in the disc and the magnetic field created by the current in the 
center electrode. Thus, the magnetic pressures created by the superposed 
fields force the charged particles of the disc to travel toward the muzzle 
of the coaxial accelerator. 

It is assumed that some of the energy acquired from the capacitors 
by the charged particles is exchanged with the neutral particles causing the 
neutral particles to diffuse into the accelerator. Also, as the charged 
particles are being driven to the muzzle of the accelerator the neutral 
particles diffuse into the region which has been vacated. 

The diffusion of the neutral particles can be described by the following 
equation: 


dn 

J = 

<L dq. 


where J is the flux of particles in the q. direction, D is the diffusion 

Q f 

dn 


coefficient, and — is the density gradient of particles. The neutral 

dq i 

gas can diffuse in all directions but will diffuse most rapidly in the direction 
of highest density gradient. Since the highest density gradient occurs in 
the negative z direction the neutral gas will expand most rapidly in a 
direction away from the muzzle. 

The neutral gas is relatively hot and will produce some ionization 
by collision processes between neutral particles. If the density is sufficiently 


68 




high to sustain an arc the electrical discharge will spread throughout the 
volume of the coaxial accelerator. The discharge will continue until the 
gas density becomes too low to sustain the arc, or until the energy stored 
in the capacitor bank has been dissipated. 
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